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ABSTRACT 

We report on Chandra observations of the bright ultra-luminous X-ray (ULX) source in NGC 3921. 
Previous XMM-Newton observations reported in the literature show the presence of a bright ULX 
at a 0.5-10 keV luminosity of 2xl0 40 erg s _1 . Our Chandra observation finds the source at a lower 
luminosity of ~ 8 x 10 39 erg s , furthermore, we provide a Chandra position of the ULX accurate to 
0.7" at 90% confidence. The X-ray variability makes it unlikely that the high luminosity is caused by 
several separate X-ray sources. In 3 epochs of archival Hubble Space Telescope (HST) observations 
we find a candidate counterpart to the ULX. There is direct evidence for variability between the two 
epochs of WFPC2 F814W observations with the observation obtained in 2000 showing a brighter 
source. Furthermore, converting the 1994 F336W and 2000 F300W WFPC2 and the 2010 F336W 
WFC3 observations to the Johnson [/-band filter assuming a spectral type of 07 I we find evidence 
for a brightening of the [/-band light in 2000. Using the higher resolution WFC3 observations the 
candidate counterpart is resolved into two sources of similar color. We discuss the nature of the ULX 
and the probable association with the optical counterpart (s). Finally, we investigate a potential new 
explanation for some (bright) ULXs as the decaying stages of flares caused by the tidal disruption of 
a star by a recoiled supermassive black hole. However, we find that there should be at most only 1 of 
such systems within z=0.08. 

Subject headings: stars: binaries — X-rays: binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are off-nuclear 
X-ray point sources in nearby galaxies wit h X-ray lu- 
minosities, L X > 1 x 10 39 - 10 42 erg s -1 (Colbert & 
|Mushotzky||1999| |Farrell et aT1|2009[ ). Their X-ray lu- 
minosities are suggestive of intermediate-mass (10 2 -10 5 
Mq) black holes (IMBHs) if they radiate isotropically 
at sub-Eddington levels as observed in stellar black holes 
and Active Galactic Nucleii (AGNs). Hence, ULXs could 
be a new kind of black hole with masses in-between the 
stellar-mass black holes found in X-ray binaries and the 
supermassive black holes (SMBHs; >lx 10 6 M ) found 
in the centers of galaxies. Thes e IMBH could be the 
building blocks of SMBHs (e.g. |Volonteri||2010[). For 
a com prehensive overview of ULXs see |Feng & Soria| 

pmj ). 

Recent studies of the luminosity functions of ULXs and 
fainter extra-galactic X-ray binaries showe d evidence for 



2011j|Mineo et al.|2012) . This, together with the identi- 
fication of an ultra-luminous spectral state under super- 
Eddington accretion rarel y observed in Galactic black 
hole X-ray binaries (e.g., jGladstone et al.||2009|), led 



to the idea that the majority ot ULXs are stellar-mass 
black holes. These stellar-mass black holes have par- 



(King et al. 


2001) anc 


sion (jBegelman 


2002) 



Nevertheless, 



a number of the 

,40, 



brightest ULXs with an X-ray luminosity > 2 x 10'" erg 

s~ x remains hard to explain as X-ray binaries. The 
higher the luminosity, the less likely it is that they can 
be explained as the high-luminosity end of the X-ray 
binary luminosity function. 

A seemingly convincing case for an IMBH is pro- 
vided ^bxihe_jwiabje j very bright ULX ESO 243-49 
X-l ' 
high 



Farrell et al. 



a break at a luminosity ~ 2 x 10 40 erg s 1 (Swartz et al. 
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2009 ). Its X-ray luminosity is too 
or a stellar-mass black hole even in the presence 
of some beaming. Given the evidence for the detection 
of a rcdshiftcd Ha emission l ine in the optical coun- 
terpart to ESO 243-49 X-l ( |Wiersema et al.| |2010"| , 
the uncertainty on its distance is reduced with respect 
to other bright ULXs. Another case for an IMBH is 
M82 X41.4+60, whereas it does not reach peak lumi- 
nosities as high as ESO 243-49 X-l, the maximum lu- 
minosity is still uncomfortably high for a st ell ar-mass 
black hole |Strohmayer & Mushotzky 2003) 

Si 



Feng & 

showed that M82 X41.4+6TTcusplays state 



Kaaretj ((2001 

changes similar to those of stellar-mass black holes when 
at sub-Eddington accretion rates. However, differences 
in interpretation exist. The temperature of the thermal 
emission is > 1 keV, which is high for an IMBH. [Okajima| 
et al. (2006) describe the spectrum as that from a slim 
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disk obtaining a black hole mass of 20-30 M Q . Similarly, 
Miyawaki et al. (2009) interprete the thermal emission 
as coming trom a thick corona and they find a black hole 
mass of < 200 M . 

A further interesting alternative explanation for some 
(bright) ULXs is that the ULX is caused by a recoil- 
ing SMBH in the f ew million years after the recoil event 
( Jonker et al.|2010 and references therein). This scenario 
is not feasible for galaxies where there is evidence for the 
presence of a central supermassive black hole from e.g. an 
AGN. The reason is that the recoiled SMBH needs to be 
replaced by a new SMBH in a subsequent merger, and 
the time for that is longer than the few million years 
lifetime of the recoiled SMBH as a ULX. 

The potential explanation of the ULXs with lumi- 
nosities above the break in the luminosity function as 
IMBHs or recoiled SMBHs warrant further investigation 
of these sources. So far, less than a dozen sources with 
L x > 10 



41 erg ; 
Wolte: 
nd IStro 



galax y, [Wolter et al.||2006 
2009 an 



1 have been found (i.e. i n the Cartwheel 
M 82 X-l, |Feng fc Kaaret 



trohmayer fc Mushotzky||20031 ESQ 24 3-49, 



Farrell et al.||2009| in JNGC 5775, |Li et al. 
in CXO J 122518.6+144545, IJonker et al.||20 



2008; and 
OfTEight 



more bright off-nuclear X-ray sources, of which 2 have 



Lx> 10 41 erg s x , were recently reported by Sutton et al 



(2012), who suggest that their properties are consistent 
with accreting IMBHs in the hard state (cf. Walton et al. 
20111. 



Here, we present Chandra, archival Very Large Array 
(VLA), and archival HST/WFPC2 and HST/WFC3 ob- 
ser vations of the br ight ULX in NGC 3921 first found 
by |Nolan et al] ( p004| called NGC 3921 X-2. NGC 3921 
is classified as a protoelliptical that was formed during 
a major merger taking place approximately 0.7 Gyr ago, 
where one of th e merging galaxies was gas rich and th e 
other gas poor (Schweizer 1996 Schweizer et al. 1996). 
Given that NG(J 3921 has an AGJN the re coiling black 
hole scenario desribed in Jonker et al. (2010) is not appli- 
cable for the ULX under study in this galaxy. We adopt 
O m = 0.3, Q A = 0.7, and H = 70 km s _1 Mpc^ 1 to 
convert the redshift of NGC 3921 to a distance measure- 
ment. 

2. OBSERVATIONS, ANALYSIS AND RESULTS 
2.1. Chandra X-ray observation 
We have obtained a 6.0 ks-long observation with the 



Chandra X-ray observatory ( Weisskopf et al.|2002[) cov- 
ering the X MM-Newton error circle of the bright ULX in 
NGC 3921 ([Nolan et al.|2004[ ) on August 14, 2011 (MJD 
55787.89362)! The Chandra observation has been per- 
formed using the S3 CCD of the Advanced CC D Imaging 



Spectrometer (ACIS) detector (|Garmire]JT997| ACIS-S). 
The observation identification (ID) number tor the data 
presented here is 13296. We reprocessed and analyzed 
the data using the CIAO 4.3 software developed by the 
Chandra X-ray Center and employing CALDB version 
4.3. The data telemetry mode was set to very faint, 
which allows for a thorough rejection of events caused 
by cosmic rays. 

Since, by design, the XMM-Newton source position 
falls near the optical axis of the telescope, the size of the 
point spread function is smaller than the ACIS pixel size. 
Therefore, we follow the method of Li et al. ( 2004 ) im- 



plemented in the CIAO 4.3 tool acis_process_events 
to improve the image quality of the ACIS data. 
Using WAVDETECT we have ident ified the source re - 



lated to the AGN in NGC 3921 (cf. [Nolan et al.|[2004| 
Using our accurate radio position of the AGIN (see below) 
we apply a boresight correction to the Chandra obser- 
vation of size 0.155" to the Right Ascencion (RA) and 
-0.079" to the Declination (Dec). For this we used the 
WCSUPDATE tool inside CIAO. We checked the position 
of the AGN after the boresight correction and it matches 
the radio position to within 0.02". However, as the ra- 
dio position is accurate to ~ 0.1" we take the latter as 
the residual 1— a boresight uncertainty. As the boresight 
correction involves the position of one source only we can 
not account for effects of rotation. 

In the boresight corrected X-ray image we detect an 8- 
count X-ray source in the XMM-Newton error circle at 
a position of RA = 11:51:07.723 and Dec=+55:04:14.98. 
In decimal degrees the position, stating in brackets the 
formal 1 sigma wavdetect uncertainty in localising the 
source on the ACIS S3 CCD, is (RA,Dec)=177.78218(4), 
+55.07083(3). The 90 per cent uncertainty in the source 
position is (RA,Dec) = (0.3",0.24"), with a total uncer- 
tainty of 0.4". This is slightly lower than the uncertainty 
of 0.5" ( 90 per cent confid ence) for low count sources 
found by Hong et al. (2005), which could potentially be 
due to the use of the sub-pixel event repositioning that 
can currently be applied. For the best position we, how- 
ever, conservatively use an uncertainty of 0.5" for the 
uncertainty in localizing the source on the CCD. Given 
the possibility that the residual boresight correction is 
due to a systematic effect we assume, conservatively, that 
the overall error on the Chandra position of the ULX in 
NGC 3921 is determined by adding the localization un- 
certainty and the boresight uncertainty (0.2" at 90 per 
cent confidence). Overall, the Chandra 90 per cent un- 
certainty on the source position is 0.7". The source is 
detected 0.3' off-axis on the ACIS S3 CCD. 

In order to verify that the source we detect in the 
XMM-Newton error circle is not due to residuals from a 
cosmic ray event, we investigated the arrival times of the 
X-ray photons from the source. Indeed, the arrival times 
of the X-ray photons do not cluster in time as one would 
expect for residual emission due to a cosmic ray e vent. 

We also d etect 4 counts from the source that INolanl 
et al.| ( |2004[ ) list as NGC 3921 X-3. The best Chandra 
p^sItioiTTs (RA,Dec) = (ll:51:03.53,+55:04:50.2) which 
corresponds to 177.76472(5), 55.08061(6) in decimal de- 
grees, where the digit in between brackets denotes the 
1 a wavdetect uncertainty in localising the source on 
the CCD. We do not detect NGC 3921 X-4 that was 

Possibly, 



also considered a ULX by Nolan et al 
the source is variable as well and it was 



(20041. 
Sefow our detec- 



tion limit. Alternatively, the detection of NGC 3921 X-4 
was spurious resulting from the incomplete modeling of 



the X MM-Newton point spread function (cf. Read et al 
201 1\ . It wa s excluded from the catalogue from Walton 



et al. 



"We 



(2( 

he 



2011 ) for that reason 

ave selected a circular region of 6 pixel («3") 
radius centered on the source position of NGC 3921 X- 
2 to extract the source counts in the energy range of 
0.3-8 keV. We limited the radius to exclude as much as 
possible the central parts of the galaxy, where extended 
emission either due to hot gas or due to unresolved faint 
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point sources may be present. Similarly, we have used 
a circular region with a radius of 80 pixels away from 
any source but on the same S3 ACIS CCD to extract 
background counts. The net number of background sub- 
tracted source counts is 7.7. From the number of back- 
ground events and the background area we find that the 
predicted number of background source photons is less 
than 0.3. For a model of a power la w with index 1.7 




cm ' s A . However, given the low number of detected 
counts, the uncertainty on this value is about a factor of 
2. 

Using standard cosmological parameters the redshift 
converts to a distance of 80 Mpc for NGC 3921, which 
makes the X-ray luminosity in the range 0.3-8 keV 
Lx ~ 8 x 10 39 erg s _1 . For comparison the 0.5-10 keV 
luminosity is the same when assuming a powerlaw spec- 
tral model with an index of 1.7 and the Galactic inter- 
stellar extinction given above. 



2.2. HST WFPC2 observations 

We have analysed archival Hubble Space Telescope 
(HST) data obtained with the Wide Field Planetary 
Camera 2 (WFPC2) in the F336W, the F555W and the 
F814W filters corresponding roughly to the Johnson U-, 
V- and Cousins /-bands, respectively. We converted the 
WFPC2 vegamag to the Johnson V- and Cousins 
/-band magnitudes using the SYNPHOT package running 
under STSDAS in iraf assuming the spectral energy dis- 
tribution of the 07 I star BZ 64. We use this procedure 
and spectral type for all filter transformations after this. 
The HST data were obtained on June 28, 1994 (MJD 
49623) and the proposal identi fication number i s 5416 . 

"~ ( fT996| 



Schweizer et al. 



Using these HST observations 

provide the V- and /-band magnitude (see their table 
3 source number 129) of a source consistent with our 
Chandra X-ray position of NGC 3921 X-2. As we will 
see below, the V and /-ba nd magnitudes o f U=2 4.03 
and /=24.02 determined by |Schweizer et aL (19961 are 
fully consistent with the values we derive below using the 
same HST data. We summed the reduced drizzled im- 
ages provided by the Space Telescope Science Institute 
giving total exposure times of 2x600 s, 2x1200 s plus 
30 s, and 2x900 s plus 30 s for the U-, V- and /-bands, 
respectively. 

We performed astrometry within midas, comparing 
the positions of stars in the field of view of the WFPC2 
against entries from the seventh release of the Sloan 
Digital Sky Surv ey (SDSS DR7) r' catalogue ( |Adelman 



McCarthy|2009 ) . An astrometric solution was computed 
by fitting ior the reference point position, the scale and 
the position angle, considering all the sources that are 
not saturated and appear stellar and unblended. We ob- 
tain a solution with 0.12" root-mean-square (rms) resid- 
uals from 15 stars well distributed on the image. The 
astrometric accuracy of SDSS DR7 is better than 75 
milli-arcseconds rms per coordinate, plus 20-30 milli- 



arcseconds due to systematic errors (Pier et al. 2003 1 . 
For the accuracy on our stellar positions we sum m 
quadrature the residuals of the astrometry with the rms 
accuracy of the catalogue and add linearly the system- 



atic uncertainty on the SDSS positions. The resulting 
positional accuracy with respect to the International Ce- 
lestial Reference System (ICRS) is 0.17" on both RA 
and Dec. at 1 a. The combined Chandra and HST 
WFPC2 positional accuracy at 90% confidence is there- 
fore 0.85" (0.7" comes from tying the Chandra astromet- 
ric frame to the ICRS, see above, and an additional un- 
certainty of 0.48" from tying the HST frame to the ICRS, 
adding both in quadrature). 

There is a bright point source present inside the 90% 
confidence astrometric error circle in all three bands 
at a position RA=11:51:07.65 and Dec=+55:04:15.2 
([RA,Dec]=177.78187(5), +55.07089(5) in decimal de- 
grees, where the value in the between brackets denotes 
the l-er uncertainty; see Fig. [I]). We used sky values 
determined local to the source of interest as the sky 
background varies substantially due to a variable amount 
of unresolved stars of the galaxy. We performed point- 
spread function fitting for the p hotometry using the HST- 
phot package (Dolphin 2000). The candidate coun- 
terpart present in the summed images gives F336W= 
23.27±0.15, F555W= 23.92±0.08, F814W= 24.01±0.11 
(U = 23.7, V = 23.9, / = 24.0). The error on the re- 
ported photometry is mainly due to spread in aperture 
correction. As the magnitudes in the U-, V—, and I— 
band depend on the spectral type assumed when using 
SYNPHOT, we do not provide errors for their magnitudes. 

Using hstphot we placed artificial point sources of 
varying magnitude in the region of 150x150 pixels 
around the Chandra position to determine the limiting 
magnitude of each of the three combined images sepa- 
rately. We detect point sources at 5-er with F336W< 
24.6,' F555W< 26.8, F814W< 25.7 (U < 25.0, V < 26.8, 
/ < 25.7). The bright optical source is not resolved in the 
WFPC2 images although the sharpness parameter from 
hstphot is at the extreme end of the range for point 
sources pointing at a nearly resolved source. 

A second set of HST WFC2 observations was obtained 
on August 2, 2000 using the F300W filter (a wide [/-band 
filter) and the F814W filter for total exposure times of 
1900 s and 260 s, respectively. Using hstphot we find 
that the F300W magnitude of the source is 22.3±0.1, and 
the F814W magnitude is 23.5 ± 0.1. In both cases the 
limiting magnitude of the exposures is ps 24. See Table [l] 
for the conversion of these F300W and F814W magni- 
tudes to Cousins / and Johnson [/-band magnitudes. 

2.3. HST WFC3 observations 

In addition to the WFPC2 HST observations, we have 
analysed archival HST data obtained with the Wide Field 
Camera 3 (WFC3) in the UVIS F336W and the F438W 
filters. The WFC3 UVIS detector has a pixel scale of 
0.04" per pixel, which is a factor of 2.5 better than that 
of WFPC2, and 0.13" per pixel in the infrared (IR) band. 

The HST WFC3 data were obtained on Oct. 3, 2010 
(MJD 55472) and the proposal identification number is 
11691. We summed the reduced drizzled images provided 
by the Space Telescope Science Institute giving total ex- 
posure times of 2300 s and 950 s for the UVIS F336W 
and the F438W filters, respectively. Furthermore, we 
analysed the observations obtained in the near-infrared 
F110W- (1565 s in total) and F160W-bands exposures 
(2518 s in total). 

We used a procedure similar to the one we used for 
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Fig. 1. — All the images are made using the 1994 WFPC2 observations. The small circle with 0.85" radius present in three sub-plots 
represents the 90% confidence astrometric error circle on the position of the ULX. North is up and East is to the left in all panels. Top left 
panel: Overview of the field of NGC 3921 in the F555W filter. The box approximates the size of the small zoomed images. Top right panel: 
The combined, 1200 s-long, HST WFPC2 F336W image (£/— band), showing the blue source near the edge of the small circle. Bottom left 
panel: Zoom-in of the combined 2430 s WFPC2 F555W (V-band) image around the position of the bright ULX in NGC 3921. Bottom 
right panel: Zoom-in of the combined 1830 s WFPC2 F814W (/-band) image around the position of the bright ULX in NGC 3921. 
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Fig. 2. — North is up and East is to the left in both panels. Left panel: Zoom— in of the combined 2300 s WFC3 F336W image around 
the position of the bright ULX in NGC 3921. The circle of 0.76" radius represents the astrometric 90% confidence region of the position of 
the ULX. Right panel: Zoom-in of the combined, 950 s— long, HST WFC3 F438W image, showing that the source at the edge of the 0.76" 
radius 90% confidence astrometric error circle is resolved into two sources. 
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the WFPC2 images to obtain an accurate astrometric 
solution for the WFC3 images. We considered all the 
sources that are not saturated and appear stellar and 
unblended. Using 11 stars on the F336W-band UVIS 
image and 7 stars on the near-infrared F110W image for 
the astrometry the fit for the reference point position, 
the scale and the position angle yielded rms residuals 
of 0.05" for the UVIS F336W-band image and 0.1" for 
the F110W near-infrared image. The astrometric so- 
lution is against the SDSS catalogue for the F336W- 
band UVIS image and against 2MASS for the near- 
infrared image. The bright point source present inside 
the Chandra error circle in the WFPC2 images is re- 
solved into two sources in the UVIS images with posi- 
tions for star 1: RA=11:51:07.647 and Dcc=+55:04:15.26 
and for star 2: RA=11:51:07.661 and Dcc=+55:04:15.36 
(star 1: [RA,Dec]=177.781862(14), +55.070906(14); star 
2: [RA,Dec]=177.781921(14), +55.070933(14) in decimal 
degrees, where the value in the between brackets denotes 
the l-er uncertainty; see Fig. [2]). The source is unre- 
solved in the near-infrared images. The combined 90% 
astrometric uncertainty of tying the HST/WFC3 images 
to the ICRS using the SDSS images and of localizing the 
two stars on the WFC3 optical frames is 0.3". The to- 
tal Chandra and HST/WFC3 combined 90% astrometric 
uncertainty is therefore 0.76". 

In deriving the photometry of the sources we used 
sky values determined local to the source of interest 
as the sky background varies substantially due to the 
presence of unresolved stars from the galaxy. We per- 
formed point-spread function fitting for the photometry 
using the DOLPHOT package. The WFC3 magnitudes in 
the VEGAMAG system for the resolved candidate optical 
counterpart are given in Table [T] The F336W magni- 
tude of the two sources combined is 23.29. A WFC3 
F336W magnitude of 23.29 gives a magnitude of 23.30 
for colors similar to those of the 07 I star BZ 64 in the 
WFPC2 F336W filter. This is fully consistent with the 
magnitude of 23.27+0.15 measured with WFPC2 in its 
F336W-band filter more than 16 year prior. The com- 
bined WFC3 F438W band magnitude of the two sources 
is 23.43. Converting the WFC3 band magnitudes to 
Johnson U— and -B-band magnitudes again using the 
BZ 64 07 I in SYNPHOT we get U x = 24.3 and B x = 23.9 
and U 2 = 24.7 and B 2 = 24.4 for the F336W and F438W 
magnitudes of the two sources separately, respectively 
and U — 23.7 and B = 23.4 for the two sources com- 
bined. 

In order to estimate the limiting magnitude of the 
WFC3 images we checked the magnitudes of the faintest 
sources detected at 5-<r, and we find approximate lim- 
iting magnitudes of F336W< 27.2 and F438W< 26.3 
(£/<27.6 and B<26.3). 

3. VLA RADIO OBSERVATIONS 

The field of NGC 3921 was observed in the X-band at 
a central frequency of 8.4 GHz with the Very Large Array 
(VLA) in A-configuration on two occasions in 2003. On 
June 19, 2003 the source was observed under project code 
AN114. The source 3C 48 was used as the primary cali- 
brator and the secondary calibrator was J 1146+539. In 
total the observation provides 14 minutes on target. On 
July 14/15, 2003 the source was observed under project 
code AN115. The source 3C 286 was used as the primary 



calibrator and the secondary calibrator was J 1148+529. 
In total the observation provides 25 minutes on target. 

Data reduction was carried out according to stan- 
dard procedures within the Astronomical Image Process- 
ing System AIPS. The field was not bright enough for 
self-calibration. When combined the two observations 
gave a root-mean-square (rms) value of the noise of 29 
/xJy/beam, and no detection at the position of the ULX 
down to a 3cr upper limit to the radio flux of 87 /iJy. This 
yields a 3er upper limit on the radio luminosity at the po- 
sition of the ULX of 5.6 x 10 36 erg s _1 for a distance of 
80 Mpc. 

The only significantly detected source corresponds 
to the center of the galaxy, at RA=11:51:06.8724(2) 
Dec=55:04:43.307(2) (J2000), at 1.8 mJy. The errors on 
those positions are purely statistical. The VLA astrom- 
etry is thought to be good to about 0.1" in reasonable 
weather conditions when in A-configuration. Therefore, 
we have used 0.1" as the fiducial \—a uncertainty on the 
VLA coordinates of the AGN in NGC 3921. 

4. DISCUSSION 

We report on a brief Chandra ob servation of the field 
of the bright ULX in NGC 3921 ( |Nolan et aLl[2004| . 
Compared with the discovery XMM-Newton observation 
~9.5 years earlier, the source 0.5-10 keV X-ray luminos- 
ity has faded by a factor of ~2.5. This virtually proves 
that the high 0.5-10 keV X-ray luminosity of 2xl0 40 erg 
s _1 detected by XMM-Newton is due to a single bright 
source and not due to the chance super position of a 
group of sources. If the latter were the case they, coin- 
cidentally, would have decreased their X-ray luminosity 
by the tim e of the Chandra observation (cf. Maccarone 
etaL|[2007 |. 

Using the accurate position we obtained using Chan- 
dra we investigate archival HST WFPC2 images in the 
F300W, F336W, F555W and F814W (2x) bands (which 
we converted to the Johnson U— , U-, V-, and Cousins 
/-bands, respectively) and archival WFC3 images in the 
F336W and F438W bands (which we converted to the 
Johnson U— and i?-band magnitudes). There is a single 
bright, blue optical source in the astrometric error circle 
in the WFPC2 images. In the WFC3 images this source 
is resolved into two sources that, when combined, have 
the same F336W magnitude in the WFC3 and WFPC2 
observations of 1994 and 2010 (but see below). With a 
F336W-F438W=-0.14, U-B = 0.3, and U-V = -0.2 
color, this source is blue (the color of source 1 and source 
2 is approximately equally blue). Furthermore, if at the 
distance of NGC 3921, with a distance modulus DM of 
34.5 excluding extinction effects, its combined absolute 
U-band magnitude, My, is —10.6. The Galactic extinc- 
tion in the direction of NGC 3921 is found to be 1.1 x 10 20 



cm (Dickey &: Lockman 19901. Using the relation of 



Fredehi & Schmitt ( |1995 ["between Nh and Ay we obtain 
Ay = 0.06. Using SYNPHOT we calculate the F336W- 
F438W color the 07 I star BZ 64 would have correcting 
for extinction. The observed F336W-F438W color is too 
blue to be explained by the 07 I star even if there was 
no extinction. The color of an 05 V star (BZ 77) may 
be consistent with the observed F336W— F438W color if 
the extinctio n is slightly lower than th e E(B-V) of 0.18 
derived from Dickey & Lockman ( 1990 ). So, we conclude 
that the source (s) making up source 1 and 2 have a com- 
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TABLE 1 

Results of the point-spread function photometry on the candidate optical counterpart(s) in the astrometric error 
circle using two HST WFPC2 (top part) and WFC3 (bottom part) observations. The Johnson and Cousins magnitudes have 

BEEN CALCULATED USING SYNPHOT ASSUMING A 07 I SPECTRAL TYPE (THE STAR BZ 64). 



Observing 


Source 


Time on 


HST Filter magnitude 


Johnson (UBV) /Cousins (I) 


date 


# 


source (s) 




Magnitude 


Jun. 28, 1994 


1+2 


1200 


F336W= 23.27+0.15 


U = 23.7 


Jun. 28, 1994 


1+2 


2430 


F555W= 23.92+0.08 


V = 23.9 


Jun. 28, 1994 


1+2 


1830 


F814W = 24.01+0.11 


I = 24.0 


Aug. 8, 2000 


1+2 


1900 


F300W= 22.3+0.1 


U = 22.9 


Aug. 8, 2000 


1+2 


260 


F814W = 23.5+0.1 


J = 23.5 


Oct. 3, 2010 


1 


2300 


F336W = 23.85+0.03 


U = 24.3 


Oct. 3, 2010 


2 


2300 


F336W = 24.27+0.04 


U = 24.7 


Oct. 3, 2010 


1+2 


2300 


F336W = 23.29 


U = 23.7 


Oct. 3, 2010 


1 


950 


F438W = 23.97+0.05 


B = 23.9 


Oct. 3, 2010 


2 


950 


F438W = 24.45+0.07 


B = 24.4 


Oct. 3, 2010 


1+2 


950 


F438W = 23.43 


B = 23.4 



bined spectrum similar to that of an early type star of 
approximately type 05. 

A direct comparison of the two WFPC2 observations 
employing the F814W filter obtained with just over 6 
years in between provides evidence for variability: in 
1994 the F814W magnitude was 24.01+0.11, whereas 
the observation in 2000 found F814W= 23.5 ± 0.1. The 



XMM-Newton observation of Nolan et al. ( 2004 ) was ob- 
tained on April 27, 2002 when the source was at its high- 
est X-ray luminosity detected so far. The optical obser- 
vation closest in time to the XMM-Newton observation 
also finds the source at a slightly enhanced luminosity, 
suggesting that one of the two sources is indeed asso- 
ciated with the X-ray source. Furthermore, comparing 
the Johnson [/-band magnitudes calculated assuming an 
07 I spectral type, the source is clearly brighter in 2000 
than in 1994 and 2010, with [/-band magnitudes of 23.7, 
22.9 and 23.7 for the observation in 1994, 2000, and 2010, 
respectively, re-enforcing the evidence for variability. 

The optical counterpart is class ified as a stellar associ- 
ation or large globular cluster by Schweizer et al. ( 1996 \ 



using the 1994 HST WFPC2 observations as we men 
tioned in Section 2.2 although they did not consider the 
[/-band observations. In the WFC3 observations the 
full-width at half maximum (FWHM) of source 1 and 
2 is sa 2.4 pixels. This yields a FWHM of <0.1", which 
at the distance of 80 Mpc of NGC 3921 is about 40 pc for 
each source. This is on the small side of the size distribu- 



the associations in IC 1613 have sizes of 40 pc (Garcia 



tion of OB associations (cf . iBresolin et al. 119981) al though 

of 



et al.|2010 ). The U—V color is also consist ent with that 



ot yoimg~T 10-50 Myr old) OB associations (|Schulz et al. 



20021 Anders k Fritze-v. Alvenslebenl 120031). in a few 



nearby ULXs, e.g. NGC 1313 X-2 and Hoimberg IX X- 



Grise et al. 


2008 


Grise et al. 


2011 



to the ULX and the other is an OB cluster the separa- 
tion between them is similar to that for M 82 X-l and 



NGC 7479 found by Voss et al. (2011) 



Optical counterparts to ULXs with luminosities below 



1 x 10 erg s are typically fainter than the absolute 
magnitude for NGC 3921 X-2, My = - 10.6. The bright- 
est have M v = -9 ( [Roberts et al.|2008b . The ULXs with 
higher luminosities are associated with brighter optical 



counterparts. For i nstance, the brig ht ULX in ESO 243 
49 has Mr = —10 ( |Soria et al.|2010 1 and th e bright ULX 
in CX O 122518.6+144545 has M g , = -10.1 (I Jonker et al. 
20101) . In the following we discuss the possible nature ot 



this bright ULX and the associated potential counter- 
part. 

4.1. NGC 3921 X-2: a ULX with a stellar-mass black 

hole? 

The two candidate optical counterparts to 
NGC 3921 X-2 are equally blue in the WFC3 ob- 
servation. Since the accretion disk color can b e similarly 
blue (e.g. van Paradijs & McClintock 1995), perhaps 
the observed light of one of the sources is dominated by 
the accretion disk. If so, the other source is dominated 
by a cluster of stars of a combined color similar to that 
of the accretion disk. The observed X-ray and optical 
variability and in particular the large variability in the 
[/-band does suggest that one of the two optical sources 
identified in the WFC3 images is indeed related to the 
counterpart to the ULX. If one of the two sources is 
due to a single object, it could also be a Lum inous Blue 
Variable (LBV; Humphreys & Davidson 1994[ ). If there 
is evidence that this source is an LBV it is prudent to 
assume (Occam's razor) that it will be the mass donor 
to the accreting black hole in the ULX. The reason is 
that both LBVs and ULXs are so rare in a galaxy that 
the probability to find a ULX and an LBV very close 
together while not physically related is the product of 
the two (low) probabilities of finding one such source in 
a certain area in a galaxy. 

On the other hand, one may be inclined to assume that 
evidence for an LBV for one of the two optical sources 
implies that the other source must be the ULX, reasoning 
that the number of LBVs in a binary orbit with a black 
hole must be rarer still than single LBVs. However, most 
O st ars are formed in a binary with another (massive) 
star ( Sana et al.|2012 ), making the scenario of finding an 
LBV with a black hole plausible. 

Finally, a scenario where both blue optical sources de- 
tected in the WFC3 images are young OB clusters is also 
valid as long as one of them hosts the ULX. When the 
ULX is X-ray bright it has to have an optical counter- 
part of magnitude similar to that of the OB association, 
explaining the variability. Therefore, the scenario where 
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the ULX in NGC 3921 X-2 is explained as a stellar-mass 
black hole accreting from a donor in a wide system such 
that the absolute magnitude of the donor plus disk rivals 
that of an OB association is entirely consistent with the 
data. 

4.2. NGC 3921 X-2: a type Iln supernova? 

As the combined i7-band magnitude from the WFC3 
source 1 and 2 measured in 2010 is consistent with being 
the same as the WFPC2 [/-band magnitude measured 
in 1994 while there is evidence for moderate brighten- 
ing in 2000, a scenario explaining the ULX as a type Iln 
super nova (cf. |Gal-Yam &: Leonard 2009 ; Jonker et al. 
20101 is strongly constrained. The scenario that still is 
(theoretically) possible but that we deem unlikely is as 
follows: the progenitor star of the supernova was not 
detected, therefore it should reside in a cluster that is 
still present. The supernova went off some time between 
1994 and 2000 and, remarkably, went undetected in the 
optical. This includes a non-detection in archival ground- 
based observations from the Jacobus Kaptein Telescope 
in 1996. The supernova has faded in the optical before 
the HST observation in 2000 nearly down to the pre- 
explosion level. The X-rays come from the supernova 
and they are now fading. 

4.3. A new scenario for the nature of some ULXs. 
Tidal disruption of a star by a recoiled black hole ? 

A galaxy like our Milky Way is predicted to have hun- 
dreds of recoiled IMBHs/SMBH s with associated star 
clusters (O'Leary & Loeb 2012). For the majority of 
systems which are expected to have had recoil veloci- 
ties <1000 km s _1 , the rate of tidal disruption events 
remains constant for about 1 Gyr after the recoil event. 
Given the expec ted rate one could have a bout a 1000 
events in 1 Gyr (|Komossa fc Merritt||2008[ ). Tidal dis- 
ruption of stars from such a bound (nuclear) cluster or 
from unbound stars of the galaxy will produce periods of 
enhanced accr etion providing X-ray luminosities above 



10 40 erg s" 1 (Komossa fc Merritt 2008 Stone fc Loeb 



2012|. The short duration Fddington-limited hares may 
well be missed by current scanning X-ray All Sky Mon- 
itors, whereas the flares may be t oo soft to be detected 
by the Swift BAT instrument (cf. Lodato & Rossi 2011 
for the predicted lightcurves and the disk temperatures). 
Theoretically, one expects the timesc ale of the decay 
of ma ss accretion rate to go as t 



-5/3 



( Roes 



1988 



Rees 



19901. Recent Chandra observations ot three potential 



tidal disruption events are consistent with this theoreti 
cal trend and an overall lifetime of approxima tely 5 years 
at L x > 10 40 erg s _1 ( |Halpern etal] |2004| ) . In order 
to determine the importance ot tidal disruption events 
from recoiled, off-nuclear black holes for the population 
of ULXs we need to estimate the number of such black 
holes that recoiled in approximately the last Gyr. The 
rate of tidal disruption events drops off steeply for longer 
time s as the loss-cone of the re coiled black hole is emp- 
ti ed ([Komossa fc Mer ritt||2008[ ) 



Patton et al. 1 2000 ) tound that over the last 1 Gyr 
about 1% ot galaxies with an absolute £?-band magni- 
tude —21 < Mb < —18 went through a merger (cf. their 
formula 32). This was derived evaluating the number of 
interacting galaxies found in the Second Southern Sky 



Redshift Survey. As the authors explain, this fraction 
is a lower limit as only pairs of relatively bright galax- 
ies have been included. Including fainter galaxies would 
increase the fraction as more interacting brigh t + faint 
galaxy pai r s wou ld be found. More recently, |De Prp7 



pris et al. (2007) found a merger fraction of 2% using 
close pair fractions and asymmetries from the Millenium 
Galaxy catalogue. 

Assuming that the black hole merger time is signif- 
icantl y less than 1 G yr such as in gas-rich mergers 



(e.g. Mayer et al.||2007 [) will provide a lower limit on the 
number ot recently (within the last Gyr) recoiled black 
holes. In the case that this merger time is much longer 
(of order of several Gyrs) the number of recently recoiled 
black holes is larger as th e galaxy merger fraction was 
high at larger redshift (e.g. |Patton et al.|200"0 Lin et al. 



20041 

Tn order to derive the number of off-nuclear recoiled 
black holes that recoiled in approximately the last Gyr 
we need to know the fraction of merging black holes that 
receive a significant kick. This fraction is debated at 
present. The magnitude of the recoil kick depends on 
the mass ratio of the two black holes before their merger, 
on the orientation of their spin axes with respect to the 
binary orbital plane and on the amount of spin. The 
presence of massive gas discs around the binary black 
hole such as potentially in gas-rich mergers can align the 
spin axes of the two black holes prior to their merger, re- 
ducing the amplitude of the recoil kick imparted durin g 
the merger ( |Bogdanovic et ah 2007 Dotti et al. 2009). 
However, fragmentation of the accretion disc around the 
binary SMBH could reduce the rate of alignment as frag- 
ments could result in short r andomly oriented accretion 
events (King & Pringlc 20071. Recent calculations seem 
to suggest that the recoil kicks might be limited to ve- 
locities lower than ss 200 km s _1 due to a grav itational 
wave "anti- kick" during the ringdown phase (Le Tiec 
et al. 2010). Similarly, under certain conditions rela- 



tivistic precession will align or anti-align the bl ack hole 
spin vectors reducing the median kick velocities ( Kesden 



et al. 2010). Recoil velocities of a few hundred km s" 
imply that many recoiling SMBHs will be retained by 
the host galaxy. Too small recoil velocities will mean 
that the recoiled black hole will sink back to the nucleus 
of the galaxy and X-ray emission due to decaying tidal 
disruption events will not be identified as ULXs. For 
now we assume that 50% of the recoiled black holes will 
have recoil velocities such that they are bound to the 
galaxy and that they appear as off-nuclear sources for 1 
Gyr. We note however, that this number depends on the 
distribution of recoil kick velocities which is uncertain. 

Putting all these estimates together we calculate the 
probability of finding a ULX caused by tidal disruption 
of a star by a recoiled SMBH per galaxy. There is 2% 
(the merger fraction) x 50% (right recoil velocity) x 1 x 
10~ 6 /yr (the tidal disruption rate per year) x 5 (the 
life-time in years above a luminosity of lxl0 40 erg s _1 ) 
leading to about 5 in 10 8 galaxies that potentially has 
a bright off-nuclear X-ray source due to the fading of 
emission caused by a tidal disruption event. 

The 1 Gyr look-back time corresponds to a redshift 
of zsa 0.08. This provides a co- moving volume of 0.16 
Gpc -3 . Given the galaxy density ~ 1 x 10~ 2 Mpc -3 
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(cf. Schechter 1976 Loveday et al. 1992 Driver et al 
2005 tor h=0.7) there are about 1.6 million galaxies in 
this volume. Therefore, there will be between 0.1-1 X- 
ray source with L> 10 40 erg s _1 due to tidally disrupted 
stars from off-nuclear, recoiled, black holes within a red- 
shift of 0.08. Note that a luminosity of L> 10 40 erg s" 1 
at a rcdshift of z=0.08 corresponds to a flux of only 
R3 7 x 10~ 16 erg cm" 2 s" 1 making it difficult to detect 
these sources in the existing shallow surveys. If the es- 
timates above are reasonable there might be at most 
one of these among the (bright) ULXs such as the ULX 
NGC 3921 X-2. 



4.3.1. Application to NGC 3921 X-2 

The decay in X-ray luminosity of a factor of 2.5 in 
NGC 3921 X-2 that we find between the XMM-Newton 
and Chandra observations ~ 3400 days apart is (too) 
small compared to the predicted drop-off from the soft 



disk X-ray s pectral component predicted by jLodato fc 
Rossil (I2011I). Even if in some proposed tidal disrup 



tion events the X-ray spectrum is found to be signif- 
icantly harder than expected, the luminosity is a bout 
a factor 200 lower 3165 days after the soft peak (Ko- 
mossa et al. 2004[ ). Given the presence of nuclear ac- 
tivity in NGC 3921 a new central SMBH must have re- 
placed the recoiled black hole, for instance following a 
recent merger, or the recoiled black hole was formed by 
a three black hole interaction event. All these constraints 
of course make the scenario of a tidal disruption event in 
a recoiled SMBH for NGC 3921 X-2 unlikely. 

4.4. Alternative explanations for NGC 3921 X-2? 

Alternative explanations for the bright ULX in 
NGC 3921 such as a background AGN or a foreground 
star or (quiescent) low-mass X-ray binary (LMXB) are 
possible in principle. For NGC 3921 X-2 to be an active 
LMXB the distance should be «1 Mpc for a luminosity 
of lxl0 36 erg s _1 typical for an LMXB. This would give 
an absolute U-band magni tude of around -1 which is in 
line w ith that of LMXBs (van Paradijs & McClintock 
19951, however, finding an LMXB in interstellar space is 
mTIikely. NGC 3921 X-2 can also be a quiescent LMXB 
at a distance of a few kpc. The absolute magnitude of the 
optical counterpart would then be pointing to an M-type 
companion, which would imply that the blue candidate 
optical counterparts have to be unrelated to the quies- 
cent LMXB. 



Schweizer et al. ( 1996 1 estimate the number of back- 
ground galaxies in the WFPC2 field of view as 1.2 x 10 -3 
per arcsec -2 . The presence of the background galaxy 



cluster Abell 1400 is taken into account in this estimate. 
The optical counterpart does coincide spatially with one 
of the tidal tails that clearly belongs to NGC 3921. As 
the (young) stellar associations in NGC 3921 trace these 
tails very well the probability of finding a background 
AGN among these stellar associations is low. 

In the case that the optical source is unrelated to the X- 
ray source the observed variability in the F814W (and U— 
band) WFPC2 observations is hard to explain. Neverthe- 
less, assuming the variability is spurious or explained in 
another way, we can provide a limit on the absolute mag- 
nitude of the counterpart of My > —7.3, M B > —8.2, 
My > -7.7 and Mj > -8.8. These limits on M a , M B , 
My and Mj have been derived using the upper limits on 
the apparent magnitude of < 27.2, < 26.3, < 26.8 and 
< 25.7, respectively and given the distance of 80 Mpc 
towards NGC 3921. 

Therefore, the source could be a ULX with a stellar 
counterpart outside an association. The upper limit to an 
optical source in this scenario translates to a lower limit 
on the ratio of the X-ray flux to the optical flux of 12. 
Whereas most CTiandra-selected AGN have a ratio of X- 



ray to optical i?-band flux of lower than 10 ( Laird et al 
20091, our limit does not rule out a background AGIN 



scenario in the case that the optical source is unrelated to 
the X-ray source. The most secure way to differentiate 
between the discussed scenarios is to obtain an optical 
spectrum of the potential counterparts. 
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